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ABSTRACT 

We present medium resolution near-infrared (NIR) spectra, covering 1.1 to 
3.4 pm, of the normal Type la supernova (SN la) SN 2014J in M82 obtained 
with the FLITEGAM instrument aboard SOFIA approximately 17 — 25 days 
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after maximum B light. Our 2.8 — 3.4 /rm spectra may be the hrst ~ 3 /rm 
spectra of a SN la ever published. The spectra spanning the 1.5 — 2.7 /im range 
are characterized by a strong emission feature at ~ 1.77 /im with a full width 
at half maximum of ~ 11,000 — 13,000 km s“^. We compare the observed 


FLITECAM spectra to the recent non-LTE delayed detonation models of Dessart 


et al. (2014) and hnd that the models agree with the spectra remarkably well in 


the 1.5 — 2.7 /im wavelength range. Based on this comparison we identify the 
~ 1.77 /xm emission peak as a blend of permitted lines of Co II. Other features 
seen in the 2.0 — 2.5 /rm spectra are also identihed as emission from permitted 
transitions of Co II. However, the models are not as successful at reproducing 
the spectra in the 1.1 — 1.4 /im range or between 2.8 /rm and 3.4 /im. These 
observations demonstrate the promise of SOFIA by allowing access to wavelength 
regions inaccessible from the ground, and serve to draw attention to the usefulness 
of the regions between the standard ground-based NIR passbands for constraining 
SN models. 


Subject headings: supernovae: individual (SN 2014J) 


Introduction 


It is widely accepted that Type la supernovae (SNe la) originate from runaway ther¬ 
monuclear burning of a carbon-oxygen white dwarf that approaches or exceeds the Chan¬ 
drasekhar mass limit. However, the evolutionary scenario (e.g., accretion from a binary 


companion in the so-called single-degenerate model Hoyle & Fowler 1960 Whelan & Iben 


1973 or merger with another white dwarf in a double-degenerate model Webbink 1984 


Iben fc Tutukov||1984| ) that gives rise to the SN and the details of the explosion mechanisms 
remain uncertain. Given the importance of SN la to distance determinations and cosmology 


see e.g., Goobar & Leibundgut 2011), it is crucial to understand the physics behind these 


explosions. One of the best ways to do this is to carry out detailed comparisons between 
high quality observational spectra of bright SNe la with physical models. Near-infrared 
(NIR) spectroscopic observations of SNe la are particularly important in this regard as they 
are not signihcantly affected by the presence of foreground or circumstellar dust that can 
seriously impact and complicate the interpretation of ultraviolet and optical spectra. In 
addition, the line blanketing opacity that dominates the UV and optical regions is lower in 
the NIR, making the interpretation and attribution of features easier. In order to reduce 
uncertainties in the model parameters and improve the accuracy of the models themselves, 
a continuous spectrum of a bright nearby SN la is necessary against which the models can 
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be compared. This includes the regions between the standard ground-based NIR pass-bands 
that are inaccessible from most observatories due to strong atmospheric absorption. 


SN 2014J in M82 was discovered on 2014 Jan. 21 by Fossey et ah (2014). Zheng et ah 


(2014) estimated the time of the explosion (‘hrst light’) to be about 7 days prior to discovery, 


Jan. 14.75 UT, while Marion et al. (2014) derived the time of maximum B brightness to be 


Feb. 01.74 UT. One of the nearest and brightest SNe la in modern history, it has been studied 
across a large range of wavelengths since its discovery, allowing the determination of a number 
of its physical parameters. Ultraviolet, optical, NIR, and mid-infrared observations have 
revealed that SN 2014J was a spectroscopically normal SN la, nearly identical to SN 2011fe, 


~ 2.0; 

Goobar et al. 

2014; 

Amanullah 

Zheng et ah 

2014 

) and relatively higher 


velocity {y > 11,000 km s absorption features in its optical and NIR spectra (Goobar 


et al. 

2014| Marion et al. 

2014). The analysis of the optical and NIR spectra presented by 

Marion et al. 

(2014 

) reveals radial stratihcation of elements produced during the explosion. 


indicative of relatively little mixing and consistent with that expected by delayed detonation 
models. Telesco et al. (2014) were able to reproduce successfully the 8 — 13 pm spectra of 


SN 2014J, which are dominated by strong lines of [Co II] and [Co III], with a model from 


Hoflich et al. (2002) for a spherical delayed detonation explosion on a Chandrasekhar-mass 


white dwarf. Based on hts to the observed gamma-ray line fluxes, Churazov et al. (2014) 


and Diehl et al. (2015) estimated that ~ 0.5 — 0.6 Mq of ®®Ni was produced in the explosion. 


in agreement with the model used by Telesco et al. (2014). 


Because SN 2014J occurred in a nearby galaxy {D 3.4 Mpc; [Dalcanton et al.]|2009[ ) 
and therefore was relatively bright, it constitutes a nearly ideal object on which to test SN la 
spectral models. Fortuitously, it also provided an opportunity to test the FLITECAM instru¬ 


ment (McLean et al. 2006) on the NASA Stratospheric Observatory for Infrared Astronomy 


(SOFIA), which can be used to hll in some of the wavelength gaps in NIR spectra of SNe la 
obtained from ground-based observations. NIR observations obtained with SOFIA are much 
less susceptible to the deleterious effects of atmospheric water vapor, whose strong absorp¬ 
tions serve to dehne the standard ground-based J (~ 1.1 — 1.35 pm), H (~ 1.5 — 1.75 pm), 
and K (~ 2.0 — 2.35 pm) bands. At altitudes > 39,000 feet (~ 12 km), typical of SOFIA 
flights, the precipitable water vapor overburden is generally between a few to tens of mi¬ 
crons whereas ground-based water vapor overburdens are typically on the order of a few 


millimeters, even at the highest and driest sites (e.g., Giovanelli et al. 2001 Radford et al. 


2008). A description of SOFIA and its instrument suite can be found in the papers by Young 


et al. (2012), Gehrz et al. (2009), and Becklin et al. (2007). In this paper we present the 


hrst spectra obtained by FLITECAM during its commissioning hights on board SOFIA. The 
NIR spectra of SN 2014J obtained with FLITECAM/SOFIA and discussed here comprise a 
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complement to the data sets presented by Marion et al. (2014) and Friesen et al. (2014). 


Observations and Data Reduction 


SN 2014J was observed with FLITECAM on four separate SOFIA flights (flight nos. 
146, 147, 148, and 149) on 2014 Feb 19, 21, 25, and 27 UT, approximately 17-25 days after 


the time of maximum B brightness, or approximately 36-44 days after Erst light (Zheng et al. 


2014). The flights originated in Palmdale, CA and the primary goal was the commissioning 


of the FLITECAM instrument on the SOFIA telescope. FLITECAM is a NIR imager and 
grism spectrograph covering the 1 — 5 micron range (McLean et al. 2006 Smith & McLean 


2008 Logsdon et al.||2014 ) and both images (see Fig. 1) and spectra of the SN were obtained. 
The observations were acquired at four different altitudes between 38,000 and 43,000 feet and 
a range of air masses. The full-width at half maximum (FWHM) of the images and spectra 
on these flights was ~ 3". Dithered images were obtained in the FLITECAM J, H, and K 
hltersj^ FLITECAM was co-mounted with the HIPO instrument during these observations, a 
conhguration that precluded observations at wavelengths longwards of ~ 4 /im, and reduced 
the sensitivity longwards of ~ 2 /rm, due to high background levels resulting from the warm 
dichroic and transfer optics. The images were reduced in a manner typical for NIR images. 
Flat helds and background sky frames were constructed from dithered observations obtained 
6' away from the target. Individual frames were flat-helded, sky-subtracted, and then shifted 
and median combined to produce a hnal image in each hlter. The frames in each hlter were 
then aligned and combined. Photometric calibrations of the images were determined from 
measurements of the star BD-l-70 587, which is present on the FLITECAM images of SN 
2014J and for which 2MASS magnitudes are available. The resulting SN photometry is 


presented in Table 1. The agreement with the values given by Amanullah et al. (2014), 


Marion et al.| ( pMl) ), and [Foley et al.| ( |2014[ ) for times near those of our SOFIA observations 
reasonably goodj^ 


IS 


The spectra were obtained by nodding the target between two different positions along 
the slit, referred to as A and B respectively. Due to the fast rotation of the held and the 
constraints of pointing the telescope, the spectra were acquired via several sequences of AB 


^The FLITECAM filter passbands are available at 
http: //www.Sofia.usra.edu/Science/ObserversHandbook/FLITECAM. html 

^ We note that, because the spectrum of SN 2014J is unlike that for any standard star, the effective 
wavelengths of the filters and therefore the color corrections could be substantial, particularly in the J and 
H bands. However, in order to facilitate comparisons with magnitudes of the SN given in the literature, we 
have not made any such corrections to the values given in Table 1. 
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observations. The slit is 60" long and the A and B positions are separated by 30". The 
low-resolntion wide slit was nsed, which resnlted in a spectral resolving power R of ~ 1300. 
The log of the spectroscopic observations is given in Table 2. The data were rednced nsing 


Redux, the SOFIA facility pipeline (Clarke et ah 2014) incorporating the fspextool software 


package, which is modihcation of the Spextool package (Cashing et ah 2004) developed for 


the SpeX NIR spectrograph at the NASA Infrared Telescope Facility (Rayner et al. 2003). 


fspextool performs non-linearity correction, AB pair snbtraction, sonrce profile constrnc- 
tion, extraction and backgronnd apertnre dehnition, optimal extraction, and wavelength 
calibration for FLITECAM grism data. Wavelength calibration was performed nsing the 
OH emission lines in the spectra of the backgronnd sky. 

In order to acconnt for the differences in the detector response at the two positions (A 
and B) along the slit, we did not combine the A and B spectra nntil after the tellnric correc¬ 
tion and flnx calibration steps. The latter steps were carried ont nsing the xtellcor_general 


package (Vacca et al. 2003), and for these pnrposes we observed an AOV star dnring each 


flight. The spectra of the AOV star were rednced in exactly the same way as for the SN. 
Again, we kept the A and B beam spectra separate. We flux calibrated the A beam spectra 
of the SN with the A beam spectra of the standard and did the same for the B beam spectra. 
To remove residual telluric features resulting from the differences in observing altitude and 
air mass between the standard star and the SN, we computed telluric absorption models ap¬ 


propriate for the observing conditions, using the ATRAN code (Lord 1992), smoothed and 


binned them to the FLITECAM resolution and sampling, and generated correction curves as 
a function of wavelength from the ratios of the telluric models for the standard to those for 
the SN. These corrections were then applied to the spectra of the SN. The telluric-corrected 
and flux-calibrated A and B beam spectra of the SN were then combined. The spectra from 
each grism setting for a given date were then merged to produce a final spectrum spanning 
the full wavelength range available. We then scaled the spectra so that synthetic photome¬ 
try yielded values that matched the imaging photometry values provided by both our own 


images and those of Foley et al. (2014), Amanullah et al. (2014), and Marion et al. (2014) 


interpolated to the dates of our observations. 


3. Results 

The extracted FLITECAM spectra of SN 2014J are shown in Fig. 2. Wavelength regions 
that are inaccessible, or difficult to observe, from the ground due to strong telluric absorption 
are marked. In Fig. 3 we present a comparison between the spectra from Feb. 19 and Feb. 
27 over the 1.5 — 2.7 /im wavelength range. The signal-to-noise ratio (S/N) for these spectra 
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ranges from ~ 10 to more than 60 per pixel at the highest flnx levels; over most of the 
wavelength range, the S/N is > 20 per pixel. 

The most striking aspects of the spectra shown in Figs. 2 and 3 are the strong emission 
featnres at 1.55 — 1.65 /rm and ~ 1.77 /xm. The former emission featnre can clearly be 


seen in the later-time gronnd-based spectra of SN 2014J presented by Marion et al. (2014). 


Althongh the latter featnre is present in gronnd-based spectra of SN 2014J (Marion et al. 


2014), and can also be seen in some of the spectra of other SN la presented by Gall et al. 


(2012), Hsiao et al. (2013), and Marion et al. (2009, 2003), the FLITECAM spectra are 


the only data that span the entire featnre prohle; the gronnd-based spectra cover only a 
portion of this wavelength region dne to the poor and variable atmospheric transmission 
between ~ 1.8 and ~ 2.0 microns. The shapes and peak wavelengths of both featnres vary 
snbstantially between the two sets of observations, which were obtained only ~ 8 days apart. 
A Ganssian fit to the 1.77 /im featnre in the Feb. 19 spectrnm yields a centroid wavelength 
of 1.768 /im and a fnll width at half maximnm of 0.080 yum (~ 13,500 km s ^); on Feb. 
27, the featnre is centered at 1.784 fim with a width of ~ 0.062 fim (~ 10,500 km s“^). 
Uncertainties are on the order of 0.003 /im (~ 500 km s“^). If the observed wavelength shift 
were to reflect a change in velocity of a specific emission featnre, this wonld correspond to 
a velocity redshift of abont 2,700 km s“^, with a decrease in the velocity width of a similar 
amonnt. Althongh the gronnd-based data presented by Marion et al. (2014) conhrm the 


shift seen in the FLITEGAM spectra, the fnll wavelength coverage provided by FLITEGAM 
in the 1.8 — 2.0 /rm region allows the change in the profile to be identihed easily. A shift 
in the featnre wavelength is not readily apparent in the NIR spectra of other SNe la shown 


by Marion et al. (2009) bnt can be discerned in the spectra of the normal SN la SN2011fe 


presented by Hsiao et al. (2013). 


Based on the model of Wheeler et al. (1998), Marion et al. (2003, 2009) identify the 


emission in the 1.7 /im region as arising from a blend of Fe gronp lines along with absorption 
from Go H. However, the strength and symmetric Ganssian shape of the ~ 1.77 fim featnre 


argne against this interpretation and in favor of (primarily) a single emission featnre. Dessart 


et al. (2014) identify this featnre, as well as other strong featnres seen in the K band spectra 


of SNe la, as emission from permitted lines of Go H. In fact, the strengthening of the Go 
lines is responsible for the well-known secondary peak in the broadband NIR photometric 
light cnrves of SNe la in their model. Their DDGIO A4D1 (hereafter DDGIO) model spectra 
exhibit a strong Go H emission featnre at ~ 1.75 /im, arising from a blend of lines, primarily 
those at 1.7463, 1.7770 and 1.8069 /im, starting at abont 25 days after explosion and shifting 
redward with time. The redshift is the resnlt of the variation in line strengths of the varions 


components. Similarly, Gall et al. (2012) identify both the 1.55 /rm and the 1.77 /im featnre 


seen in the NIR spectra of the normal SN la SN 2005cf and three other SNe la at 9-12 days 
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past maximum B band emission, as permitted lines of Co II, based on comparisons with the 


W7 model of Nomoto et ah (1984) 


In Fig. 3 we compare our combined 1.5 — 2.7 /im spectra of SN 2014J at 18 and 26 days 
past B maximum (36 and 44 days past explosion) with the DDCIO model spectra]^ computed 
for dates close to those of our FLITECAM observations. The DDCIO A4D1 model is a set 
of one-dimensional (i.e., spherically symmetric), time-dependent, non-LTE calculations for 
a delayed detonation explosion of a Chandrasekhar-mass white dwarf at a distance of 10 
pc. The model includes non-thermal effects, non-local energy deposition, “huge” model 
atoms for iron group elements, and forbidden lines of both intermediate mass elements 


and iron group elements; descriptions of the model can be found in Dessart et al.| (2014), 


Blondin et ah (2013), and references therein. We shifted the model spectra to account for the 
redshift of M82 (n = 203 km s“^). We also reddened the model spectra to acconnt for both 
foregronnd reddening, due to both the Milky Way and M82, and circumstellar reddening 
and scattering, using the extinction values and the power law reddening/scattering curve 
advocated by Amanullah et ah (2014) (see also Foley et ar]|2014 ). We then scaled the model 
spectra to match the overall flux levels of the observed data; a single scale factor for each 
model spectrum was derived from a least-squares £t to the observed spectra. The strongest 


emission lines of Co II and [Co III], as predicted by the model (Blondin et ah 2015), are 
identified near the bottom of each plot. (The wavelengths of NIR Co II lines can be found 
in the tables given by Marion et ah 2009 and Gall et ah 2012 ) As can be seen, the model 


spectra match the observed spectra extremely well in this wavelength regime, although the 
time past explosion for the best-htting model to the earlier of the two observed spectra is 


somewhat discrepant. For an explosion date of 2014 Jan. 14.75 UT (Zheng et ah 2014), the 
FLITECAM spectra obtained on Feb. 19 correspond to ~ 36 days after hrst light, and the 
spectra obtained on Feb. 27 correspond to ~ 44 days past hrst light. The model spectrnm 
that matches the earlier data set best corresponds to 33.15 days past explosion. For the later 
data set, the model spectrum corresponding to 40.11 days past explosion provides the formal 
best £t to the data; however, the model spectra for 44.12 and 48.53 days past explosion £t 
the observed data in the region around the 1.77 pm feature considerably better. 


Since the model spectra are computed for a distance of 10 pc, the scale factors needed 
to match the model spectra to the observed hux levels yield estimates of the distance to 
SN2014J and therefore provide additional constraints on the models, and the predicted 
absolnte huxes, that are independent of the overall spectral shape. Farther, the distances 
derived from scaling the model spectra to observed spectra on different dates and over 


^ Kindly provided to us by L. Dessart and S. Blondin 



























different wavelength ranges provide a test of the self-consistency of the models. We hnd 
that the distances derived from the scale factors are well within the range of values for M82 
{D = 3.2 — 5.5 Mpc with most values between 3.2 — 3.9 Mpc; for example Karachentsev & 


Kashibadze 2006 Dalcanton et ah 2009 Tutui & Sofue 1997). The best £t reddened model 


spectra shown in Fig. 3 yield distances of 3.7 Mpc (for the model corresponding to 33.15 
days past explosion) and 4.0 Mpc (for the model corresponding to 44.12 days past explosion); 
the model spectra for 40.11 and 48.53 days past explosion yield distances of 4.2 Mpc and 
3.5 Mpc, respectively. The agreement of the values, both with each other and independent 
distance estimates, is fairly good, especially in light of the (~ 5 — 10%) uncertainties in 
the absolute flux levels of the observed spectra resulting from the scaling needed to match 
the interpolated photometric values. Finally, the comparison with the models also serves to 
illustrate the relative insensitivity of NIR observations to the effects of reddening, even for 
an object such as SN 2014J for which Ay ~ 2.0 mag (|Amanullah et ah 2014 Foley et ah 


2014). 


The remarkably good agreement between the observed spectra and the DDCIO model 
spectra indicates that this model captures the basic elements of the physics involved in the 
formation and evolution of the observed NIR spectra of SNe la. The comparison between 
the model spectra and the observed spectra in the H and K band further supports the 


conclusions of Gall et ah (2012) and Friesen et ah (2014), who suggested that the post¬ 
maximum spectra of normal SNe la, including that of SN 2014J, longward of ~ 1.5 /rm 
are dominated by broad emission from permitted lines of Co II, but calls into question the 
previous identification of Co II absorption features at high velocity in the K band spectra of 


other SNe la by Marion et ah (2009). It also suggests that the identihcation of the 1.77 /xm 


feature as emission from Co II, dominated by the 1.7770 /xm transition, is probably correct, 
despite the fact that the best fitting model spectrum predicts a somewhat wider and shifted 
line compared to the observations on Feb. 19. The discrepancy between the observed and 
predicted line shape suggests that the relative strengths of the Co II lines contributing to 
the 1.77 /xm feature are incorrect in the model. The 1.7463 /xm line is far stronger in the 
model than the observations on this date indicate. As demonstrated in Fig. 4, htting the 
emission feature with multiple Gaussians, constrained to have the same velocity width and 
centroids fixed at 1.746, 1.777 and 1.807 /xm, reveals that the observed line shape in the 
data from Feb. 19 can be successfully reproduced with relative line intensities in the ratio of 
~ 0.45 : 1 : 0. To reproduce the line shape observed in the Feb. 27 data requires of ratios of 
relative intensities of ~ 0.25 : 1 : 0.65. Therefore, the observed shift in the 1.77 /xm feature 
is due to the weakening of the 1.7463 /xm line and the strengthening of the 1.8069 /xm line. 
The line widths also decrease from 
and the optical depth decreases. 


12, 000 km s ^ to ~ 7, 000 km s ^ as the ejecta expand 
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Following the suggestion of Gall et al. (2012), we used the width of the unblended Co II 


lines at 2.36 /im and 2.46 /im to derive an estimate for the extent in velocity space of the Co¬ 
rich core. The weaker Co II line at 2.36 /rm seen in the Feb. 27 data yields a velocity extent 
of ~ 7,000 km s“^, in agreement with the line widths estimated above for the components 
of the 1.77 /im feature. The stronger Co II 2.46 /im line, the full width of which is generally 
inaccessible from the ground, yields a velocity extent of the Co core of ~ 11,000 km s“^ 
for both observation dates (Feb. 19 and Feb. 27), in good agreement with the location of 
the iron-group zone in models such as W7 and DDCIO. Since almost all of the Co results 
from the decay of ^®Ni, the velocity extent of the Co should reflect the distribution of ^®Ni 
immediately after the explosion. As can be seen in Fig. 1 of |Dessart et al.| (2014), the DDCIO 
model predicts that 99% of the ®®Ni should be located within 11, 600 km s ^ ( |Blondin~et 
al. 2013, 2015). Our estimated velocity range is in good agreement with this value, as well 
as with the values found by Gall et al. (2012) for SN 2005cf and SN 2002bo at ~ 30 — 40 
days past maximum B light (slightly later than our observations of SN 2014J at ~ 18 — 26 
days past maximum B light). 

Despite the success of the DDCIO model in reproducing the H and K band spectra of 
SN 2014J, the model does not match as well at other wavelengths. In Fig. 5 we show the 
comparison between the FLITECAM spectra obtained on Feb. 25 and 27 and the DDCIO 
model over the entire 1 — 3.5 /rm range. Clear discrepancies between the data and the model 
are readily apparent in the J band, the region around 2 /rm, and longwards of ~ 3 /im. The 
disagreement between the model spectra and the observed spectrum in the 2 /rm region at 44 
days past explosion is somewhat surprising given how well the model matches the spectrum 


in this wavelength region earlier in the temporal evolution. Although Friesen et al. (2014) 


suggest that the 2 /rm feature seen in their spectra of a sample of SNe la (including SN 
2014J) is due to emission from [Ni II] at 1.939 /xm, the DDCIO model indicates that this 
feature is due primarily to [Co III] at 2.002 pm. In fact, very little ®®Ni should remain at 
this time, as 44 days past explosion is ~ 7 times longer than the 6 day half-life of ®®Ni, 


and therefore Friesen et al. (2014) speculated that the [Ni II] emission might arise from the 
stable isotope ^®Ni. However, the observed emission feature exhibits a peak at ~ 1.99 pm, 
a wavelength which is far longer than that expected if [Ni II] were the primary component 


(~ 1.94 pm). Friesen et ah (2014) noted a similar wavelength discrepancy for this feature 
(relative to the predicted emission from [Ni II]) but attributed the wavelength difference to 
a redshift of the observed emission due to the combined effects of an asymmetrical explosion 
and Doppler shifts. Attributing the emission to [Co III] instead seems to be a far more 
satisfying explanation, especially given the success that spherical models have in reproducing 
the spectra. Nevertheless, the DDCIO model predicts this [Co III] feature decreases in 
strength far faster than the observations indicate; the emission peak is still very strong 44 
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days past explosion, whereas the models predict the feature should have faded completely 
by about day 40. 


In Fig. 6 we present a comparison between the DDCIO model spectra and our J band 
spectra. While the model seems to agree reasonably well with our observed J band spectra 
at the earlier time in the evolution (~ 36 days past explosion), the mismatch between the 
observations and the model at later times (~ 44 days past explosion) is particularly bad. 
The most prominent characteristic of this wavelength region is a relatively strong emission 
peak at ~ 1.29 pm, present in both of our J band spectra (Fig. 1). The model spectra, 
particularly for late times in the temporal evolution of the SN, do not adequately reproduce 
this feature. Similar emission peaks can be seen in some of the late time spectra presented 
by Marion et al. (2009) and Gall et al. (2012), who attribute the emission to a blend of Fe II 
lines. Again, as for the 1.77 pm and 2.46 pm features, the significantly higher atmospheric 
transmission afforded by SOFIA observations allows this 1.29 pm feature to be discerned 
and measured easily in the FLITECAM spectra. The comparison between the observed 
2.8 — 3.4 pm spectrum and the model spectra is shown in Fig. 7, where a large discrepancy 
is also apparent, despite the low signal-to-noise ratio of the data. (The low signal-to-noise of 
this spectrum is partially due to the additional background from the warm transfer optics in 
the dual instrument conhguration.) The FLITECAM data indicate that broad emission is 
present between ~ 3.0 pm and 3.4 pm, which can provide additional constraints on models 
of SNe la. This spectral region is very difficult to observe from the ground due to strong 
telluric absorption, and our data may represent the hrst ~ 3 pm spectra of a SN la ever 
obtained. 


Conclusions 


By providing access to wavelength regions that are difficult or impossible to observe from 
the ground, SOFIA can contribute significantly to studies of nearby bright SNe. We have 
demonstrated this with our FLITECAM observations of a normal SN la, SN 2014J in M82, 
about 36 — 44 days after explosion, which have fully revealed a strong emission feature at 
1.77 pm. This feature shifted redward in wavelength and simultaneously narrowed between 
our two sets of observations separated by ~ 8 days. Because this feature extends into a 
wavelength regime that is difficult to observe from the ground, the identification, width, and 
the shift of the feature have never before been appreciated. Comparison with the non-LTE 


delayed detonation DDCIO A4D1 model of Dessart et al. (2014) suggests that this feature is 


most likely due to a blend of permitted emission lines from Co 11. The variation in location 
and width is most likely due to changes in the relative strengths of the component lines. 
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Other peaks seen in the K band spectra are also due to permitted Co II lines emitting 
at the systemic velocity, as opposed to highly redshifted absorption features as suggested 
previously. The width of the Co II 2.46 /xm, a transition which is generally inaccessible from 
the ground, indicates that the Co-rich core extends to about 11,000 km s“^ in agreement 
with predictions from models. The Dessart et al. (2014) model matches the observed spectra 
remarkably well across the 1.5 — 2.6 pm range, and scaling the model spectra to match the 
observed flux levels yields distances consistent with previous estimates for M82. However, 
the model provides a much poorer match to the observed spectra in the J band, and at 
2 pm and ~ 3 pm around 44 days after explosion (~ 26 days after B maximum). The J 
band spectra are characterized by a strong broad emission peak at ~ 1.29 pm, which is not 
reproduced in the model, particularly at late times in the temporal evolution . Similarly, 
the observed spectrum between 3.0 pm and 3.4 pm, the hrst SN la spectrum ever obtained 
in this wavelength range, exhibits broad emission that is not seen in the model. The region 
around 2 pm exhibits a strong emission feature of [Co III] that decreases in strength far less 
quickly than the model predicts. Our results serve to draw attention to the usefulness of 
the NIR regime, including the regions between the standard ground-based NIR passbands, 
for constraining SN models, due in part to the relative insensitivity of the NIR to the 
effects of reddening. The wavelength coverage afforded by FLITECAM on SOFIA allows 
the identihcation and study of the evolution of features not previously recognized in ground- 
based spectra of SNe la. 
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Fig. 1.— Three color image of the SN 2014J held obtained with FLITECAM. Blue is J, 
green is H, and red is K. The arrow identihes SN 2014J, which is readily apparent as the 
point source SW of the nucleus and in the plane of the disk of M82. North is up and East 
is left. 
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Fig. 2.— FLITECAM flux-calibrated spectra of SN 2014J on the four different observing 
dates. Regions of strong telluric absorption from the ground (transmission < 20%) are shown 
in dark gray, while regions of moderate tellnric absorption from the gronnd (transmission 
< 80%) are shown in light gray (cf., for example, Rayner et ah 2009). 
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Fig. 3.— Comparison between the observed H + K band spectrum of SN 2014J obtained 
with FLITECAM (in gray) on 2014 Feb 19 UT (top) and 2014 Feb. 27 UT (bottom) and 
the CMFGEN model DDCIO A4D1 of Dessart et ah (2014) for three different times past 
explosion. The model spectra have been scaled to match the flux levels of the data. Dashed 
lines are model spectra that have been reddened due to the effects of dust in both the Milky 
Way and M82. Triangles denote the locations of lines of Co II and [Co III] predicted to be 
strong in the model. 
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Fig. 4.— Decomposition of the observed 1.77 /im feature in the spectra obtained on 2014 
Feb. 19 (top) and 2014 Feb. 27 (bottom) into three Gaussians hxed at the locations of the 
strongest Co II lines (1.7463 /im - blue, 1.7770 /im - green, and 1.8069 /rm - red). The solid 
black is the sum. The dashed lines mark the centroids of the observed feature. The observed 
shift in the feature is due to the changing line strengths of the three components. 
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Fig. 5.— Comparison between the observed spectrum of SN 2014J obtained with FLITE- 
CAM on 2014 Feb. 25 (2.8 — 3.4 /rm) and Feb. 27 (JHK bands) UT (in gray) and the 
CMFGEN model DDCIO A4D1 of Dessart et ah (2014) for three different times after explo¬ 
sion. The model spectra have been scaled to match the flux levels of the data. Dashed lines 
are model spectra that include reddening due to both Milky Way and M82 dust. 
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Fig. 6.— Comparison between the observed J-band spectrum of SN 2014J obtained with 
FLITECAM on 2014 Feb. 21 UT (top) and 2014 Feb. 27 UT (bottom) and the CMFGEN 
model DDCIO A4D1 of Dessart et ah (2014) for three different times past explosion. The 
model spectra have been scaled to match the flux levels of the data in this wavelength range. 
Dashed lines are model spectra that include reddening due to both Milky Way and M82 
dust. 
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Fig. 7.— Comparison between the observed 2.8 — 3.4 /im spectrum of SN 2014J obtained 
with FLITECAM on 2014 Feb. 27 UT and the CMFGEN model DDCIO A4D1 of Dessart 
et ah (2014) for three different times past explosion. The model spectra have been scaled to 
match the flux level of the data in this wavelength range. 
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Table 1. FLITECAM/SOFIA Photometric Observations of SN 2014J 


MJD 

UT Date 

Phase 

(days) 

Filter 

Mag 

56707.3 

2014 02 19 07:11 

17.5 

J 

10.97 ±0.06 

56707.3 

2014 02 19 07:19 

17.5 

H 

9.64 ±0.03 

56707.3 

2014 02 19 07:29 

17.5 

K 

9.54 ±0.04 

56709.3 

2014 02 21 07:22 

19.5 

J 

10.93 ±0.05 

56709.3 

2014 02 21 07:32 

19.5 

H 

9.56 ±0.04 

56709.3 

2014 02 21 07:43 

19.5 

K 

9.49 ± 0.04 


Note. — Phase is the number of days past the time of 
maximum B light (MJD 56689.8; Marion et al. 2014). 



Table 2. Log of FLITECAM/SOFIA Spectroscopic Observations of SN 2014J 


MJD 

UT Date 

t - to 

t — tsmax Grism Wavelengths texp 

A"exp Altitude ZA 



(days) 

(days) (microns) (s) 

(feet) 


56707.4 

2014 Feb 19 08:52 

35.6 

17.6 

H_C 

1.500 

- 1.718 

300 

4 

41,000 

41 

56707.4 

2014 Feb 19 10:26 

35.7 

17.6 

Hw_B 

1.675 

-2.053 

150 

8 

41,000 

37 

56707.5 

2014 Feb 19 10:50 

35.7 

17.7 

Kw_C 

1.910 

- 2.276 

150 

8 

41,000 

41 

56707.5 

2014 Feb 19 11:28 

35.7 

17.7 

Kl^ 

2.270 

- 2.722 

180 

2 

41,000 

43 

56707.5 

2014 Feb 19 11:35 

35.7 

17.7 

K1_A 

2.270 

- 2.722 

120 

2 

41,000 

43 

56709.4 

2014 Feb 21 09:35 

37.7 

19.6 

J_B 

1.141 

- 1.385 

180 

4 

43,000 

40 

56709.4 

2014 Feb 21 09:53 

37.7 

19.6 

K1_A 

2.270 

- 2.276 

120 

6 

43,000 

40 

56713.4 

2014 Feb 25 10:17 

41.7 

23.6 

LM_C 

2.779 

-3.399 

150 

4 

40,000 

35 

56713.4 

2014 Feb 25 10:31 

41.7 

23.6 

LM_C 

2.779 

-3.399 

200 

12 

40,000 

35 

56715.3 

2014 Feb 27 06:20 

43.5 

25.5 

J_B 

1.141 

- 1.385 

180 

4 

38,000 

34 

56715.3 

2014 Feb 27 06:37 

43.5 

25.5 

H_C 

1.500 

- 1.718 

180 

4 

38,000 

33 

56715.3 

2014 Feb 27 06:53 

43.5 

25.5 

Hw_B 

1.675 

-2.053 

180 

4 

38,000 

33 

56715.3 

2014 Feb 27 07:10 

43.6 

25.5 

Kl^ 

2.270 

- 2.722 

180 

4 

38,000 

33 

56715.3 

2014 Feb 27 07:27 

43.6 

25.5 

Kw_C 

1.910 

- 2.276 

180 

4 

38,000 

33 


Note. — UT date and time given are at the beginning of an exposure set. The exposure time fexp is 
the time for a single exposure, and A^^exp is the number of exposures obtained, t — to is the number of 
days past to, the hrst light date estimated by Zheng et al. (2014); t — tsmax is the number of days past 
maximum B light estimated by Marion et al. (2014). ZA is the average zenith angle of the target during 
the observations. 



